Rapid state purification protocols for a Cooper pair box 
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We propose techniques for implementing two different rapid state purification schemes, within 
the constraints present in a superconducting charge qubit system. Both schemes use a continuous 
measurement of charge (2) measurements, and seek to minimize the time required to purify the 
conditional state. Our methods are designed to make the purification process relatively insensitive 
to rotations about the x-axis, due to the Josephson tunnelling Hamiltonian. The first proposed 
method, based on the scheme of Jacobs [Phys. Rev. A 67, 030301 (R) (2003)] uses the measurement 
results to control bias (z) pulses so as to rotate the Bloch vector onto the i-axis of the Bloch sphere. 
The second proposed method, based on the scheme of Wiseman and Ralph [New J. Phys. 8, 90 
(2006)] uses a simple feedback protocol which tightly rotates the Bloch vector about an axis almost 
parallel with the measurement axis. We compare the performance of these and other techniques by 
a number of different measures. 
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I. INTRODUCTION 

Superconducting charge qubits (Cooper pair boxes) are 
a promising technology for the realisation of quantum 
computation on a large scale 0, 0|. For conventional 
fault-tolerant quantum computing, the quantum states 
should have a high level of purity, preferably being as 
close to a pure state as possible. When the qubit is 
coupled to an environment it is subject to decoherence, 
which will typically result in a completely mixed state 
0]. However, a qubit initially in a completely mixed 
state can be purified by measurement. Here we consider 
continuous measurements, which can be considered as a 
rapid succession of 'weak measurements' 0]. This gives 
rise to stochastic 'quantum trajectories' that 'unravel' 0] 
the average density operator evolution described by the 
Markovian master equation 0]. The quantum trajectory 
is for the conditional density operator, conditioned upon 
the specific measurement record that was obtained in a 
given experiment. 

In the Bloch sphere representation, pure qubit states 
lie on the surface of the sphere, with mixed states being 
in the interior, and the completely mixed state being at 
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the centre. When a weak measurement is performed on 
the qubit, the effect is to pull the state (on average) to- 
wards one of the poles on the measurement axis (which 
we will take to be the z axis for simplicity) . This pull cor- 
responds to an increase in the average purity over time. 
An infinitely fast measurement would instantly project 
the state to one of the poles, as in the quantum Zeno ef- 
fect 0]. However real measurements are never infinitely 
fast. Moreover, for a charge qubit it can be difficult to 
connect and disconnect a strongly coupled (fast) measur- 
ing device without introducing additional environmental 
noise. Thus it is necessary to consider measurements giv- 
ing a finite rate of purification. 

Since purification takes a finite time, it makes sense 
to consider whether the information in the measurement 
record can be used to change the process of purification 
via feedback. The use of such quantum feedback tech- 
niques to increase the purity of conditioned states has 
been the subject of a number of recent studies 0,0,0]. 
Jacobs showed that the maximum increase in the average 
purity occurs when the qubit Bloch vector is rotated onto 
the x-y plane (the plane perpendicular to the measure- 
ment axis) , after each incremental measurement . This 
strategy is optimal in the sense of maximizing the fidelity 
of the qubit with some fixed pure state at a given final 
time, as has recently been shown using rigorous tech- 
niques from control theory [Io| . In addition, this feed- 
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back protocol is deterministic because even though the 
conditioned density operator evolution is stochastic in 
general, the stochastic term is proportional to the pro- 
jection of the Bloch vector along the measurement axis 
Hereafter, in this paper, this protocol is referred to 
as 'ideal protocol I'. 

Although ideal protocol I performs best in the sense 
just defined, Wiseman and Ralph have recently shown 
that there are reasons to consider the (conceptually) op- 
posite approach, namely keeping the Bloch vector aligned 
with the measurement axis [9[. They show that this re- 
sults in the majority of qubits reaching a given level of 
purity earlier than in ideal protocol I. In fact this strategy 
is optimal in the sense of minimizing the expected time for 
a qubit to reach a given level of purity (or fidelity with a 
fixed pure state) [HI]. This is achieved at the expense of 
having this time be stochastic (unlike ideal protocol I). 
Specifically, the distribution of qubit purification times 
is heavily skewed, with a long tail of low purity values. 
Hereafter, in this paper, this protocol is referred to as 
'ideal protocol II'. 

This paper addresses a complication which occurs 
when one attempts to apply either of these two schemes 
to a specific model of a voltage-controlled charge qubit. 
A superconducting charge qubit consists of a supercon- 
ducting island (also known as a Cooper pair box) coupled 
to a bulk superconductor via a small capacitance and a 
Josephson weak link junction. The Josephson junction, 
which allows the tunneling of Cooper pairs onto and off 
the island, normally has limited controls. Although the 
tunnelling rate can often be modified in experimental sys- 
tems the Josephson tunnelling energy provides an 
avoided level crossing between the two qubit energy lev- 
els, and maintaining this minimum energy gap minimises 
the risk of thermal excitation of the system. Close to this 
avoided crossing the energy states are formed from super- 
positions of the quantised charge states (q = 0, 2e) that 
act as the computational basis for this qubit. The tun- 
nelling gives rise to a a x Hamiltonian corresponding to a 
rotation about the rr-axis. As the junction energy should 
not be zero, the qubit Bloch vector is continually in mo- 
tion. The effect of the Hamiltonian evolution often dom- 
inates the evolution of the system under the action of the 
continuous measurements. This means the Bloch vector 
can neither be stopped near the x-y plane or z-axis nor 
can the direction of rotation be reversed. However, the 
applied voltage bias allows some control over the z-axis 
rotations due to the a z term in the Hamiltonian. This 
voltage bias is more commonly expressed as an effective 
biasing charge n g , [TH [T3| ■ 

In this paper we study the purification of charge qubits 
(Sees. II and III) using quantum feedback, based on both 
ideal protocols discussed above. We suggest mechanisms 
that could provide near optimal purification rates in the 
presence of more realistic feedback constraints than those 
considered in the ideal protcols previously studied, 0, [§] ■ 
For the protocol based on ideal protocol I we show that 
good rates for the increase of the average purity should 



be achievable using a constant Josephson energy and 
applying controlled voltage bias field pulses to create 
a z-rotation to rotate the state vector onto the cc-axis 
(Sees. IV and V). We refer to this as practical protocol 

I. This is advantageous in two ways: firstly the x-axis 
is trivially on the x-y plane which satisfies the rapid pu- 
rification condition, and secondly the effective radius of 
the cc-rotations is reduced, so the vector remains near to 
the plane even when the control pulses are not accurately 
applied. Next, we will show that the Bloch vector can be 
constrained to the region near the measurement axis (z- 
axis), by using a strong voltage bias field applied at the 
correct moment to encourage tight radius orbits around 
an axis almost parallel to the z-axis (Sees. VI and VII). 
This decreases the average time for the qubit to purify, as 
in ideal protocol II. We refer to this as practical protocol 

II. In Section VIII we conclude with a brief summary of 
the results. 



II. SYSTEM MODEL 

A superconducting charge qubit, shown in Fig.[T] con- 
sists of a small island of superconducting material con- 
nected via a Josephson junction of tunneling frequency 
v to a bulk superconducting electrode, where v = Ic/2. 
The electrode supplies a voltage bias, which can be ex- 
pressed as an effective charge q g . The island is also ca- 
pacitively coupled to a grounded electrode to supply a 
common reference. For simplicity, we ignore the dynam- 
ics of the biasing circuitry [1 31 ] - 

The Hamiltonian of this system is 
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Here 9 is the superconducting phase difference across the 
Josephson junction expressed in units of the flux quan- 
tum $o = h/2e. There exists a commutation relation 
between the conjugate variables of charge and phase, 
[q, 8] = —ih. The capacitance C q is the effective qubit 
capacitance calculated from the three physical capaci- 
tances [l4j . Cj, C g and C p , 
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(See Appendix [Al for values). At low energies this Hamil- 
tonian can be approximated by using just two states. 
Using n g = q g /2e for the effective number of Cooper 
pairs induced by the bias voltage, the Hamiltonian in 
the charge basis is [ll|, 03 ■ 
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FIG. 1: (Color online) The simple charge qubit model is an is- 
land (orange box) coupled to a bulk bias electrode via a single 
Josephson junction whose frequency is known. The qubit is 
controlled via a bias voltage applied across the device. There 
are three major capacitances, Cj is the capacitance of the 
Josephson junction, C g couples the island to a grounded elec- 
trode, and C p is the parasitic capacitance observed between 
the bulk electrodes. 



The first term may be discarded as the identity matrix 
does not affect the dynamics of the system, but is in- 
cluded initially for comparison with Eq. |T]). The second 
(<7 2 ) term shows that the applied voltage bias field con- 
trols the rotations about the z-axis. When n g = 0.5 the 
rotations are halted; when n g > 0.5 the qubit rotates 
in one direction, and when n g < 0.5 the direction is re- 
versed. The third and final (a x ) term is caused by the 
Josephson junction. For a single Josephson junction at 
constant temperature and magnetic fields, the frequency 
of rotation around the x-axis is fixed by manufacture — 
we take the frequency to be 10GHz, in line with exper- 
imental values [l5j. This frequency is equal to the min- 
imum splitting (n g = 0.5) shown by Figure [5] and it is 
vital that a sizeable separation is maintained to preserve 
the two distinct states and suppress the effect of thermal 
fluctuations. In some implementations it is possible to 
use a flux-controlled Josephson junction to vary v fl6l ]. 
However the method proposed in this paper uses the volt- 
age bias alone to apply the feedback, thereby simplifying 
the control system. 



III. CONTINUOUS MEASUREMENT 

For scalable quantum computing [I?) it is necessary 
to work with pure states. For qubits, this means states 
on the surface of the Bloch sphere. Both pure states 
(surface of sphere) and mixed states (inside the sphere), 
can be written as a 2 x 2 density matrix p. The on- 



FIG. 2: (Color online) Energy level structure of the qubit. 
The avoided crossing caused by the 10GHz Josephson junc- 
tion creates an energy gap between the energy states. This 
separation is kept relatively large to reduce problems caused 
by thermal excitations or other extraneous noise. 



diagonal elements represent the populations of, and the 
off-diagonal elements represent the coherences between, 
the charge states. The impurity of the qubit state can 
be quantified by the following [3]: 



L = 1 - P = 1 - Tr 



(4) 



Any pure state has an impurity of zero, and the com- 
pletely mixed state has an impurity of 0.5. 

One way to increase the purity of a state is through 
measurement, that is to extract some information about 
the state of the system. Assuming continuous measure- 
ment of charge (cr z ), the conditioned state Pr of the 
charge qubit obeys a stochastic master equation [18| with 
the following form 



dp c = - r [H, p c ] alt - - [a z , [a z , p c ]] dt 



-r- (vzPc + PcOz - 2 (cr z ) p c ) dW. (5) 

The first term is the Hamiltonian evolution, with H given 
by Eq. The second term represents the back-action 
of the measurement (parametrized by 7). In the absence 
of Hamiltonian evolution, this causes a deterministic de- 
cay of the mixed state towards the z-axis. (We have 
assumed that there are no other sources of decohcrencc 
for simplicity). The third term is due to conditioning 
upon the measurement result. It is stochastic, with dW 
being a Wiener increment (3I. [l9|. That is, in every time 
interval of duration dt, dW is an independent Gaussian 
distributed random variable, with zero mean and vari- 
ance equal to dt. This stochastic term depends upon the 
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particular unravelling considered [J, [l||, which de- 
pends upon the measurement scheme. In this case we are 
conditioning the state upon a continuous 'current' which 
is different in every run of the experiment and which is 
given by 

^Tv[p c (t)a z ]dt + J^-dW(t) (6) 



I{t)dt = -±Tr[p c (t)a z ]dt + 

This is a current in the generalised sense used in quantum 
optics and other areas, such that I(t)dt is dimcnsionless. 
If this measurement result is ignored then one simply av- 
erages over the last term in the stochastic master equa- 
tion j5]). This yields the deterministic master equation 
given by the first two terms of Eq. §5§ . 

Rather than evolve equation ((5|), we implemented the 
simulations using Bloch coordinates v = [x,y,z] T which 
is equivalent to using the density matrix formalism, how- 
ever the positional coordinates are easier to visualise. 
The equations for the incremental changes in the Bloch 
coordinates due to continuous weak measurement dv = 
[dx, dy, dz] T can be found in Appendix iBl in addition the 
rotations due to the non-zero Hamiltonian acting on the 
Bloch vector must be included. The equations for dx, dy 
and dz are then numerically integrated over time. 

As the system evolves under continuous measurement, 
it tends to be pulled towards the surface of the Bloch 
sphere as information (the measurement record) is ob- 
tained. In the absence of Hamiltonian evolution this 
Bloch vector will be aligned with the measurement (z) 
axis, and the system will evolve stochastically towards 
one or other of the two poles defined by this axis through 
the Bloch sphere. If the Hamiltonian evolution is in- 
cluded and if 7 is relatively small, the Bloch vector will 
rotate under the action of the Hamiltonian and be only 
weakly perturbed by the measurement interaction. The 
information extracted by the measurement is dependent 
upon the orientation of the Bloch vector with respect to 
the measurement axis @. This means that manipulat- 
ing the Hamiltonian by external controls can affect the 
way that the purity of the system increases. It is this 
that forms the basis of the rapid purification protocols 
discussed in the following four sections. 



Tn(e) is defined implicitly by Lu{Tu) = e, where @, [1] 
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This integral can be solved numerically Q, but for long 
times (small e)an analytical approximation gives Tu ~ 
ln(e~ 1 )/47 For shorter times (larger e) Tu is bigger 

than this expression Q. In general, Tu can be used to 
define the speed up factor for a given test method, 
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Ideal Protocol I 
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The ideal protocol I [3] rotates the Bloch vector onto 
the plane orthogonal to the measurement axis to max- 
imise the increase in the average purity during each in- 
cremental step. In our situation, this means rotating 
onto the x-y plane. This protocol eliminates the stochas- 
tic contribution to the evolution of p c , as can be verified 
from the stochastic equations in Appendix [Bj Thus the 
impurity equals the average impurity, which decays ex- 
ponentially: 



U(t) 
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From this, the time taken to reach impurity e is 
T, = ln(e-V2)/8 7 , 



(9) 



(10) 



and for minimizing this time this an exceptional protocol. 
Thus the maximum speed up factor is 2, in the limit of 
very small e. 

It would be extremely difficult to apply these instant 
and perfect control fields to a practical qubit. In addi- 
tion, for our superconducting charge qubit there is also 
the continual motion of the state vector due to the non- 
zero Josephson junction energy. The protocols discussed 
below address these issues. 



B. Flux-controlled Hamiltonian Feedback 



IV. FEEDBACK PROTOCOLS I 

In this and the following section we are concerned with 
protocols for maximizing the increase in the average pu- 
rity, as in ideal protocol I. Thus we need a baseline by 
which to compare the various methods. This baseline is 
given by ideal protocol II — a situation in which the ideal 
feedback controls cancel any Hamiltonian evolution and 
the qubit Bloch vector is allowed to drift stochastically 
towards the poles. This appears to be the worst in terms 
of the time T[e) taken for the average impurity L to drop 
to a given level e. For the ideal protocol II this function 



Although this is not feedback, natural cc-axis rotations 
take the Bloch vector through the x-y plane (Fig. [3]), so 
there is still an improvement over having no Hamiltonian 
evolution at all. The spiral path only momentarily passes 
through the x-y plane so it does not experience the full 
benefit of the rapid purification protocol. However, it is 
possible to utilise a method which uses a flux-controlled 
Josephson junction to manipulate the a;-axis rotational 
frequency [16[ . The algorithm slows the qubit whilst the 
Bloch vector is close to the x-y plane and then hastens the 
passage through the z-axis poles, maximising the time 
close to the x-y plane. The control is purely via the a x 
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term (modulating the Josephson tunnelling frequency, v) 
and always maintains a significant energy gap to suppress 
the effects of thermal fluctuations. This benefit of this 
approach is significant [l6| but not as close to ideal pro- 
tocol I as the following approach, ft is important to note 
that Figure [3] should not be interpreted as the average 
position of the Bloch vector, as the stochastic measure- 
ment noise causes random initial phases for the rotation, 
and as such the average position of the Bloch vector is 
the centre of the Bloch sphere for all time. This is also 
true for Figures HJ El and [TU1 which are only provided to 
illustrate the feedback concepts. 




1 -1 

FIG. 3: (Color online) Under continual s-rotation and the in- 
fluence of weak measurement, the time evolution of the Bloch 
vector takes a spiral path. The weak measurements increases 
the purity, as the purity is a direct measure of the Bloch vec- 
tor length, the length of the vector is increasing. As there 
is only one axis of rotation, the vector length corresponds to 
the radius of the orbit. Hence the spiral path found in the 
j/z-plane only. 

(Note that this should not be interpreted as the average posi- 
tion of the Bloch vector, as the stochastic measurement noise 
causes random initial phases for the rotation, and as such the 
average position of the Bloch vector is the centre of the Bloch 
sphere for all time. This is also true for Figures 3J [S] and I10I 
which are only provided to illustrate the feedback concepts.) 



C. Practical Protocol I 

The first novel algorithm proposed in this paper at- 
tempts to use finite duration voltage bias pulses to ro- 
tate the Bloch vector repeatedly on to the x-axis, tak- 
ing a screw- like path (Fig. |3J\.) , we refer to this as prac- 
tical protocol I. The a;- axis is of particular interest as 



it is invariant under x-rotation. Therefore if the vector 
can somehow be positioned close to the x-axis, it should 
remain close to the x-y plane, even in the absence of 
further control pulses. This is why the x-axis is an at- 
tractive target in the presence of continuous x-rotation 
due to Josephson tunnelling. However, the effect of the 
weak measurement is to pull the Bloch vector towards 
the poles, so the Bloch vector will be gradually pulled 
away from the x-axis in a growing spiral path. Therefore, 
to successfully return the Bloch vector to the x-axis, a 
simple control scheme has been devised which utilises a 
finite duration Hamiltonian proportional to a z to return 
the vector to the x-axis within a half cycle. The feed- 
back process triggers when z exceeds a particular thresh- 
old 

^Limit (Fig. On triggering, the controller applies 
a bias field (z-axis rotation) of the required amplitude 
and duration. An advantage of this approach is that the 
control field does not need to be continually altered. 

After the Bloch vector has reached the x-axis, the bias 
field is removed, so the qubit only experiences the con- 
stant x-rotation. This creates a distorted spiral path 
(similar to that in Fig. [3]) which will once again expand to 
exceed ZLimit, where the feedback will trigger. The over- 
all effect is to constrain the magnitude of z to ZLimit, so 
that the Bloch vector remains relatively near the x-axis. 

Ignoring the effects of weak measurement during the 
feedback pulse, it is possible to determine the pulse am- 
plitude and duration analytically. Consider a point in 
the xz-plane with z = ZLimit (Fig. [5]) and use a 7r-rotation 
about an axis (dotted arrow) to finish exactly on the x- 
axis, by rotating from the very top to the very bottom 
of the circular path about this axis. The angle a of the 
axis can be calculated using elementary geometry from 
•^Limit-axis and the distance from the centre of the Bloch 
sphere (|a| = 1 — 2L). The latter can be obtained from 
the conditional density matrix p c (which represents the 
current state of knowledge of the system) . We find 

This angle is bounded as follows: 

cw = \ sin" 1 (1) = 45° (x = 0) (12a) 

"min = ^ Sin_1 ( Z Limit) (H = 1) (12b) 

To implement this control strategy it is necessary to 
determine the x and z angular velocity components. To 
simplify matters we assume that the Josephson junction 
angular frequency is fixed, tu x = v. Then we find we 
require 

lo z = uj x ta,na (13) 

the required bias value can then be obtained from equa- 
tion (fT5)) below. Figure [5] shows the relation between the 
measured z value and the application of feedback as a 
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FIG. 4: (Color online) (A) Pictorial representation of the pro- 
posed practical protocol I which attempts to return the Bloch 
vector to lie on the x-axis. This is achieved using pulsed rota- 
tions around the z-axis, the combined rotational frequencies 
of the x and z rotations define a plane of rotation about an 
arbitrary axis. The timing and duration is selected to move 
the vector from the top of the circular path to the bottom (x- 
axis). This quick corrective feedback is performed whenever 
the weak measurement process pulls the vector past a pre- 
determined threshold (here ZLimit = 0.5), creating a screw- 
like path. The dashed vertical lines represent the gradual 
spiralling observed in Figure [3] (B) An example trajectory 
showing the effects of the stochastic continuous measurement 
noise for the first application of the feedback pulse train. The 
noise corrupted feedback path finished above the target x- 
axis, where the second spiral is highlighted in red. 




FIG. 5: (Color online) The applied feedback creates a rotation 
about the axis (dotted arrow) defined by a, this angle can 
be calculated through a measurement, or estimate, of z and 
either the purity, P or impurity, L. The z rotation frequency 
is adjusted and pulsed to ideally take the Bloch vector to 
(|a|,0,0) within a half cycle. 

function of time for the first two applications of the feed- 
back protocol. For this example ZLimit = 0.333 and the 
other values can be found in Appendix [X] The bias pulse 
duration r, required to perform the 7r-rotation about the 
a axis is 



2 y/iV x 2 +UJ Z 2 

The relation between the z-axis rotational frequency 
and bias is a linear relationship described by the following 
equation: 

For our feedback mechanism presented we find that the 
maximum u> z is equal to the Josephson junction fre- 
quency, which means there is a limit on the size of bias 
which should be applied. This is actually a favourable 
constraint as applying a bias field substantially larger 
than n g > 0.75 increases the risk of accessing an un- 
wanted third charge state [HI, HH . Another requirement 
of the control system is being able to halt the ^-rotations 
or at least slow the rotations significantly by setting the 
n g close to 0.5. It is expected that both of these require- 
ments should be achievable. The bias control range to 
compensate for a system with a constant 10 GHz cc-axis 
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rotation and the capacitances provided in Appendix IA1 
is: 



Z Limit =0.333 




0.5000 sC n g ^ 0.5323. 



(16) 



The pulse train featured in Figure [S] shows a decreasing 
trend in the voltage bias amplitude (n g ) and an increase 
in the pulse duration (r). This increase is due to the 
slower z-axis angular velocity at the latter stages. The 
values of n g and r tend to steady state values defined 
by a m in, (Eq- (|12bp ). Of course, as Figure 2J3 shows, the 
evolution of the Bloch vector will not be as smooth or pre- 
dictable as that depicted in the illustration Figure [4j\, in 
addition the control pulse cannot be applied or removed 
instantaneously There will always be some control de- 
lay and some uncertainty as to when the Bloch vector 
is likely to exceed the threshold value. The stochastic 
terms will make the prediction of the evolution uncertain 
and, consequently, the timing of the pulses will contain 
an uncertainty. To mitigate potential problems in the 
timing of the control pulses, the evolution of the Bloch 
vector can be filtered further to reduce the effect of the 
noise (the evolution described by equation (5) already 
represents a filter of the information extracted from the 
qubit). However, numerical calculations including tim- 
ing errors show that the accuracy of the control pulses is 
not critial to the performance of the purification proto- 
col (see Figure 10 below). As long as the Bloch vector 
is reasonably close to the threshold value and the pulse 
takes the Bloch vector back to the vicinity of the x-axis, 
then the majority of the available purification increase is 
still obtained. This is good news for a practical imple- 
mentation of this protocol because it demonstrates that 
the approach is robust to errors in the control system. 
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FIG. 6: (Color online) Idealised timing diagram for the feed- 
back pulses with ZLimit = 0.333 and using the values found 
in appendix [A] The topmost graph shows the ideal evolu- 
tion of z as a function of time, with the dashed and solid 
line segments corresponding to those found in Figure[3]4.. The 
dashed lines indicates no feedback (n g = 0.5), and so the si- 
nusoids are rotations about the rr-axis due to the tunnelling 
across the Josephson junction (Fig. [3J). When ZLimit = 0.333 
is exceeded, the feedback is applied, and the Bloch vector is 
7r-rotated from the top of the tilted plane to the bottom such 
that 2 = 0. The process is repeated when z again exceeds 
ZLimit, the pulses eventually tending to the same shape set by 
fEas. ll2bll4li5l 



V. RESULTS I 

As mentioned in the previous section, finding the peak 
value of z could prove difficult under noisy conditions. 
However, we have found that finding the peak is not re- 
quired if the measurement of the system is sufficiently 
weak such that the threshold ZLimit is not far exceeded. 
Using, a simple threshold on the z value derived from p Cl 
the resulting performance is still close to the ideal. 

Figure [7] shows the improvement in the average purifi- 
cation rate, or (for convenience) the 'speed-up', [7(. This 
improvement is measured relative to the minimum aver- 
age purification protocol, the case of free measurement 
evolution with no Hamiltonian (Eq. 0), which forms 
the baseline performance. The shape of the graph indi- 
cates that the performance increase is not constant for 
all values of purity, with maximum gains at high purity 
(the final part of the time evolution) . It should be noted 
that this graph is also independent of the measurement 
strength 7. Equation §E§ defines the improvement S as 
a function of remaining impurity L, and exact analytical 
solutions this equation are non-trivial, however it is pos- 
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FIG. 7: (Color online) The graph of purity rate improvement 
shows that the practical feedback protocol proposed in this 
section (dashed line) achieves similar performance to that of 
the ideal method (solid line), which required instantaneous 
feedback. In addition, the feedback mechanism described in 
section IVI CI (dotted line) has little performance gain, as ex- 
pected. It should also be noted that using no feedback in this 
system still yields a minor improvement (dash-dot line) over 
having no Hamiltonian evolution at all 



sible to invert the run-averaged transients graphically, or 
use piecewise approximation [22j ■ 

The ideal protocol I requires ideal and instantaneous 
feedback which would be very difficult to achieve in prac- 
tice. This motivated considering the no feedback case, 
for which the Bloch vector continually rotates about the 
x-axis. The constant rotation takes the Bloch vector 
through the equatorial plane twice per cycle, this momen- 
tarily approximates the ideal feedback conditions. This 
creates a minor increase in the purification rate, the dash- 
dot line (Fig. [T]). 

When the feedback routine detailed in section llVl is ap- 
plied to the qubit with a best possible threshold value for 
^Limit; an almost optimum increase is achieved (dashed 
line) using a single practical control field. Hence we have 
shown that is is possible to create a control strategy to 
gain a significant amount of the ideal purification rate. 

The performance of the algorithm depends on the value 
of the ZLimit threshold, (Fig. [H]). We find good perfor- 
mance over a relatively large range, 0.2 < ZLimit < 0.6. 
This implies the system can tolerate inaccurate and noisy 
thresholding, with only a minor penalty in the perfor- 
mance. For low values of ZLimit the thresholding per- 
forms poorly due to the system operating under the 'noise 
floor': due to the measurement noise during the control 
pulses triggering often recurs before the Bloch vector is 
returned to the vicinity of the x axis. For large values 
of ^Limit the trigger for applying the feedback is only ac- 
tivated at a late stage, and so the speed-up is severely 



FIG. 8: (Color online) This graph of the average purification 
rate improvement as a function of ZLimit shows the dependence 
on the thresholding at three different levels of purification. 
Setting ZLimit ^ 1 is equivalent to deactivating the feedback 
as the Bloch vector is constrained inside a unit sphere, al- 
ternatively using a small value threshold (zLimit = 0.1) also 
yields lower performance due to the inherent measurement 
noise. The horizontal dashed line is the theoretical maximum 
performance increase given by Ideal Protocol I at L = 0.001. 
The flattened peak for L = 0.001 near this 'envelope' indi- 
cates a high tolerance of inaccuracies, with triggering between 
z = 0.2 and z = 0.6 being quite sufficient with an optimal 
value of ZLimit ~ 0.333 as indicated by the vertical dashed 
line. 




Q- 10 20 30 40 50 60 70 80 90 

Phase delay 



FIG. 9: (Color online) Using the optimal value for the thresh- 
olding, in this case ZLimit = 0.333. This graph of the average 
purification rate improvement as a function of the time delay 
starting from the top of the Bloch vector orbit (expressed as 
a phase delay with frequency F x ) suggests there is a reason- 
ably large window of opportunity to apply the control fields 
and still achieve near optimal results. Indeed this is a test of 
synchronisation, it could also be possible to allow the Bloch 
vector to make complete rotations and later apply the fields 
within this window. 
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reduced. 

Increasing the frequency v of the Josephson junction 
aids the feedback procedure as the increased tangential 
velocity of the rotation reduces the time for the Bloch 
vector to reach the x-axis. This is advantageous as the 
measurement noise disrupts the path taken (an example 
is provided in Figure [4)3), therefore the reduction in the 
time a point travels increases the probability of reaching 
close to the desired destination, the x-axis. In addition, 
the relatively large Josephson junction energy gap re- 
duces the possibility of thermal excitations between the 
two energy levels. 

The anticipated disadvantage of using such high fre- 
quencies would be timing problems, although the control 
pulses should be feasible as the required 7r-pulscs arc of- 
ten used in quantum information processing. If time de- 
lays are a problem, then the qubit state can be allowed 
to rotate through several complete cycles if the radial 
growth of the Bloch vector is sufficiently small. Figure [9] 
shows the effect of delaying the application of the feed- 
back as a phase angle from the top of the Bloch vector 
orbit, at the Josephson frequency of 10GHz for an op- 
timum ^Limit ~ 0.333. Ideally, the feedback should be 
applied immediately (0°) but it can be seen that there 
is an allowable delay with a small decrease in expected 
performance. 

To summarize, we have shown in this section that 
the proposed method achieves a near optimum improve- 
ment in purification, and further improvement could be 
achieved by an increase in Josephson junction frequency. 
It can be clearly observed in Figures [8] and [9] that the 
method is quite tolerant of errors, and it is found to be 
sufficient to rotate the qubit back to the vicinity (rather 
than exactly on) the x axis. 



VI. FEEDBACK PROTOCOLS II 

In contrast to ideal protocol I, the ideal protocol II 
Q maximizes the stochastic terms in equations (|B1[) by 
keeping the Bloch vector on the measurement axis. 

By using feedback to rotate the Bloch vector to the 
measurement axis, it has been shown recently that the 
majority of qubits reach a given level of purity faster 
than by using ideal protocol 10. It is the existence of 
rare but extremely long purification times which makes 
the average purification time longer for this scheme than 
that for the deterministic ideal protocol I. 

For a qubit which has either Hamiltonian evolution 
around the measurement axis or no Hamiltonian at all, 
this protocol requires no controls to implement., as is not 
a problem and the Bloch vector naturally diffuses along 
the measurement axis between the two possible outcomes 
(the poles). However, there is an issue when a Hamilto- 
nian takes the Bloch vector away from the measurement 
axis, meaning the effect of the system measurement will 
be less than ideal. Unfortunately, this is the case for our 
superconducting charge qubit, where the continual x-axis 



rotation takes the Bloch vector away from the z-axis and 
passes through the x-y plane before returning to the z- 
axis (Fig. [3]). This helps for the purpose of maximizing 
the rate of purity increase of a single qubit, but is harm- 
ful for the purpose of minimizing the average time for a 
qubit to reach a given purity. 

A. Ideal Protocol II 

Assuming ideal and instantaneous controls, it would 
be possible to apply feedback to rotate back perfectly to 
the z-axis after each measurement. If this were possible 
the purity would evolve as per equation ([7]). 

B. Flux-controlled Hamiltonian Feedback 

A tunable flux-controlled Josephson junction can be 
used to change the x-axis rotational frequency, slowing 
the passage of the Bloch vector near the poles and speed- 
ing it through the equatorial plane (l6| . This is the dual 
of the method described in section IIVB1 

C. Practical Protocol II 

It is anticipated that there would be difficulties in 
adapting the previously described method, for rotating 
to the x-axis, to rotate to the z-axis instead. Due to 
the constant x-axis rotations and the switchable z-axis 
rotations, the nature of the problem is not symmetrical. 
Whenever the z-axis rotations are removed, the Bloch 
vector will still rotate about the x-axis therefore taking 
the vector away from the required z-axis. Unless the ex- 
perimental apparatus can measure, process and apply a 
correcting control field within a fraction of a cycle, the 
application of feedback will be futile as a complete cycle 
about the x-axis will have been made anyhow. 

To solve this, we propose a very simple feedback pro- 
tocol where the Bloch vector rotates about a tilted axis 
almost parallel to the measurement axis (Fig. [TOR) . 

Initially, no z control is applied and the Bloch vector 
is allowed to rotate about the x-axis. The effect of the 
weak measurement is to pull the Bloch vector towards 
the poles, and a spiral growth results. This initial period 
allows an experimentalist to detect a sizeable peak or 
trough in the z measurement record corresponding to the 
phase of the oscillation in z, indicating when the Bloch 
vector is near to the z-axis. When at this point, if the 
threshold value (zLimit = 0.333) is exceeded, the strong z 
control is applied, creating a rotation about an axis that 
should be almost parallel to the z-axis. If the initial de- 
tection is completed successfully, the Bloch vector should 
be near the z-axis and will now travel in a tight spiral 
close to the z-axis (Fig. fTUR). Alternatively, if the Bloch 
vector was somewhere near the y-axis, for example due 
to an initial delay, the orbital path about this tilted axis 
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FIG. 10: (Color online) (A) Pictorial representation showing 
the intended path of the Bloch vector. The feedback control 
is initially off, so that the Bloch vector continually rotates 
and grows around the a;-axis, this is to allow the peak value 
of z to become more distinct. Once the Bloch vector has 
exceeded a threshold and is at a maximum or minimum, the 
high frequency z-rotations are applied which locks the Bloch 
vector to the measurement axis. (B) An example trajectory 
showing the effect of measurement noise on the path of the 
Bloch vector. Whilst the noise is significant, the overall shape 
of the stochastic trajectory is similar to the conceptual path 
shown above. The shape confirms that the Bloch vector is 
confined to the vicinity of the measurement axis (z-axis) as 
required. 




FIG. 11: (Color online) Angle a as a function of bias, n g . This 
illustrates the operating range of the two techniques described 
in this paper. The green region indicates the necessary an- 
gles and biases for practical protocol I, the range of which is 
intrinsically constrained by the maximum value of a = 45°, 
(0° < a < 45°). The solitary dashed line illustrates a possible 
bias value for generating the tilted axis of practical protocol 
II. Notice that only one value is required. We require a to 
be as close to 90° as possible, however if the bias is set too 
far from 0.5, there is a risk of accessing an unwanted third 
state. For large bias values the gains in a are actually quite 
minimal. 

would be much wider and fail to coincide with the actual 
z-axis. One could instead simply apply this <j z control 
from the start of the purification process, to ensure that 
the qubit always remains close to the z-axis. However, 
there is little benefit in applying the control at the early 
stages because the performance gain close to the centre 
of the Bloch sphere is small. 

The definition of a as the angle of the axis of rotation 
between the x and z axes, still holds for this method. 
Equation (|17|) defines a as a function of Bloch sphere 
coordinates. Here we express it in terms of the system 
frequencies: the constant Josephson junction frequency 
u> x and the bias control frequency lu z : 

-MS)- <17) 

The magnitude of the bias field should not be too large 
or the next charge state may be accessed and the two 
state approximation would be violated. Figure [TT1 shows 
how a varies as a function of the bias control, plotted 
until n g = 0.75 as this is halfway between the charge 
states. A bias value of n g = 0.70 is chosen in the simula- 
tions below to reduce the possibility of accessing a new 
state. This gives a = 82° (see Fig. [TTJ) . Increasing n g 
gives minimal gain in a as the angle asymptotes to 90°. 
Ideally we would have a = 90° but we find an angle of 
82° gives acceptable performance. 
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VII. RESULTS II 

Figure [12] shows that the performance gain for the 
practical protocol II is near that of ideal protocol II. 
This indicates the desired operation: a practical imple- 
mentation of rotating to the z-axis despite the presence 
of the constant x-axis rotations due to the qubit tun- 
nelling. Thus considering Figs. [7] and [T2l we see that 
both the practical protocols I and II, emulate the ideal 
counterparts. 



Ideal protocol I 
Practical protocol I 
No feedback 
Practical protocols 
Ideal protocol II 




10 10 
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FIG. 12: (Color online) The graph of the improvement, in 
comparison with ideal protocol I, in achieving the time of 
first passage shows that the practical protocol proposed in 
this section (blue line) achieves similar performance to that 
of the ideal protocol II (black line), which required perfect in- 
stantaneous feedback. In addition, the practical protocol de- 
scribed in section II V CI (red line) has little performance gain, 
as expected. It should also be noted that using no feedback 
in this system still yields a minor improvement (green line) 
over having no Hamiltonian evolution at all 

To demonstrate the usefulness of rotating the Bloch 
vector to the z-axis, we examine the statistical distri- 
butions for the five options described in this paper. In 
Figs. [T51 and [TTl the following notation is used: A - Per- 
fect rotations to the z-axis (Ideal protocol II 0), B 
Rotations about an axis almost parallel to the z-axis 
(Practical protocol II), C - a x tunnelling Hamiltonian 
only (No feedback), D - Rotating to the x-axis (Prac- 
tical protocol I), and finally, black dashed line - Perfect 
rotations to x-y plane (Ideal protocol 10)- The compo- 
nent values and measurement strength used to simulate 
the following histograms can be found in Appendix lAl 

Figure [T51 depicts the distributions of the times of first 
passage, that is, the time for the stochastic impurity 
to fall below a given impurity of 1 x 10 ~ 3 . Each his- 
togram comprises 10 5 independent simulation runs sep- 
arated into 100 bins. Figure [T5R confirms that for ideal 
protocol II, the majority of runs reach the required pu- 
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FIG. 13: (Color online) Plotting the distribution of times to 
reach a given impurity of 1 x 10 -3 clearly shows that the ma- 
jority of simulated runs (100,000 runs) of the other protocols 
reach the target impurity before the deterministic ideal pro- 
tocol I, (indicated by right most dotted line), which greatly 
reduces the average impurity. Due to the deterministic na- 
ture of ideal protocol I, it guarantees that all qubits reach a 
set impurity in a given time. However it has been suggested 
that the poor average performance of rotating to the z-axis 
is due to the existence of a few important and extraordinar- 
ily long duration runs (> 20ns simulation time) which have a 
significant effect on the calculation of the average. 



rity at times earlier than that predicted by ideal protocol 
I. The modal value for ideal protocol II (grey dashed line 
in all the plots) is 4.0ns as opposed to 10.4ns for ideal pro- 
tocol I. It is apparent that the poor average performance 
is due to the existence of a small number of extraordi- 
narily long duration runs which have a significant effect 
on the average purity. Note that the many runs which 
do not reach the required purity by the end of the simu- 
lation (T max = 20ns) and are included in the histogram 
marked as '> 20ns'. 

Figure I13B shows the remarkable similarity between 
the ideal and practical protocol II distributions, with the 
modal values in close alignment. This would indicate that 
rotating about a tilted z-axis is a near optimal approach 
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FIG. 14: (Color online) The histograms (50,000 runs) of the 
final impurity at the simulation end time of 7.5ns, demon- 
strates the large range of final impurities given ideal protocol 
II (A), allowing access to better purity in a shorter time at the 
risk of worse purity on some occasions. Interestingly, when 
perfect feedback is replaced with practical protocol II (B), 
the range of very small impurities (below 10 ) is not acces- 
sible. In addition, when using practical protocol I we can (for 
these example values) be confident that by 7.5ns almost all 
impurities will be at least smaller that 5 x 10 -2 , indicated by 
the arrow. 



for rotating to the z-axis, given the presence of a constant 
x or y rotation. 

When no feedback is applied, there is a constant rota- 
tion about the x-axis. This rotation momentarily passes 
the Bloch vector through the z-axis and the x-y plane, 
creating a compromise between the two ideal purification 
methods (Fig.[T3p). 

Rotating the Bloch vector to the x-axis, (practical pro- 
tocol I as described in section llV C|) . yields a distribution 
of smaller variance. Indeed it should be noted, as indi- 
cated by the arrow, no runs exceeded the maximum sim- 
ulation time of 20ns. The modal value is closer to that 
of ideal protocol I (10.4ns). The non-zero variance is due 
to the need to allow the Bloch vector to grow and rotate 
about the x-axis; the non-zero z component (Fig|?j3) con- 



tributes measurement noise to the purity [see Eq. (|B1[) ] . 

Figure Q3] shows the distribution of impurities at t = 
7.5ns, which is the time at which the impurity under ideal 
protocol I reaches 5 x 10~ 3 . Each histogram contains 
50,000 runs separated into 50 logarithmically spaced bins. 

Comparing Figs. [14"R and [HP we see a dramatic dif- 
ference in the spread of values by many orders of magni- 
tude. The deterministic natures of both ideal protocol I 
and the reduced stochasticity of the more practical proto- 
col I can easily be observed. Of particular interest is the 
area corresponding to high impurity indicated by the ar- 
row. In Figure HMD this region is mostly unoccupied, but 
the other three histograms which do not employ protocol 
I have high occupancy. This implies that although these 
three methods can potentially reach very low impurities, 
it is done at the risk of ending with a high impurity. 

Interestingly, Figure [T4B follows a similar profile to 
Figure [14"R until the impurity is of the order 10 -5 , when 
smaller impurities become inaccessible. This is due to 
a mushrooming effect which creates an end cap to the 
expected path of the Bloch vector. The end cap occurs 
whenever the Bloch vector is near a pole at the surface 
of the Bloch sphere and is due to the weak measure- 
ment noise. This can be further explained by examining 
Eqs. (|Bip for weak measurement in the Bloch sphere rep- 
resentation when the Bloch vector is near a z-axis pole. It 
can be seen that as |z| approaches one, the random con- 
tribution of the Weiner increment becomes much larger. 
As the x and y values are non-zero (due to the off-axis ro- 
tations removing the Bloch vector from the z-axis), the 
'large' random changes in x and y combined with the 
constant rotation due to Hamiltonian evolution makes 
it naturally improbable that the Bloch vector will settle 
exactly on the z-axis. Hence, in practice access to the 
smallest impurities may be difficult without increasing 
the measurement strength 7 relative to v. 



VIII. 



CONCLUSIONS 



We have considered two techniques for rapid state pu- 
rification for use with a model superconducting charge 
qubit with a single control field and continuous weak 
measurement. We show that near optimal results can 
be obtained using a realistic implementation of feedback 
control. For practical protocol I, the feedback is simple 
to calculate and uses constant amplitude 7r-pulses that 
are applied for time scales that are comparable with the 
natural period of the qubit evolution. In addition, as the 
maximum angle between the rotation axis and the x-y 
plane (a) is 45°, the range of bias controls is confined to 
a small range of values close to the default bias condition. 
For practical protocol II, the feedback control is simply 
maintained, once triggered. 

If an experimenter wished to ensure that the major- 
ity of qubits reach the same level of purity at the same 
time, ideal protocol I or the more practical implementa- 
tion of rotating to the x-axis as described in section HV CI 
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should be used (practical protocol I) . Alternatively, if the 
objective is to maximise the number of qubits attaining 
a given level of purity, the experimenter should choose 
practical protocol II described in section fVI CI and based 
on ideal protocol I. As the techniques need only keep the 
Bloch vector close to the ideal conditions, the practical 
protocols are expected to be robust to a variety of con- 
trol errors (e.g. time delays, magnitude of bias, pulse 
duration) . 

Both of the practical protocols (I and II) described 
in this paper operate with continuous rotations about 
the rr-axis, which are generated by the constant tun- 
nelling frequency of the Josephson junction (which is 
the realistic scenario for practical superconducting charge 
qubits). Numerical simulations demonstrate that both of 
the practical protocols perform well and are not adversely 
affected by this constraint on the controls allowed in the 
system. In fact, the constant a x term arising from the 
tunnelling is essential to the correct operation of prac- 
tical protocol I. If these continuous rotations are either 
naturally occurring or can be applied, it opens the pos- 
sibility of implementing these purification techniques in 
other systems that contain such Hamiltonian evolution. 



APPENDIX A: TABLE OF VALUES 

Values are constant and consistent for all simulations, in 
line with experimental values quoted in reference [l5j . 



TABLE I: Components 





Description 


Typ. 


v/2-k 


Josephson junction energy 


10GHz 


Cj 


Josephson junction capacitance 


500aF 


c B 


Qubit-Grounded Bulk capacitance 


0.5aF 


C p 


Electrodes parasitic capacitance 


l.OaF 


7 


Measurement strength constant 


75 x 10 6 



APPENDIX B: WEAK MEASUREMENT IN THE 
BLOCH SPHERE REPRESENTATION 



Here we state the Cartesian equations for the random 
incremental changes in x, y and z for each time step 
dt due to a continual weak measurement process with 
measurement strength 7. These equations can be found 
in reference p^ |. and they can be derived by working 
through equation ([5]) using the Pauli matrix identities 
and equating the resulting density matrix elements with 



the Bloch vector coordinate equations: x 
y = Tr {p c (T y } and z = Tr {p c o\J. 



dx 
dy 
dz 



-(Ajdt + z^/8jdW)x 
-(4<ydt + z-y/8rfdW)y 

(1 - Z 2 )y^dW 



L 



2 2 2\ 

x -y - z ) 



Tr {pccr x }, 



(Bla) 
(Bib) 
(Blc) 

(Bid) 



This set of simultaneous stochastic differential equations 
is not trivial to solve. Indeed, equation ([7]) is the special 
case where x = y = 0, and yet yields an integral that 
appears to have no analytical solution Q ■ 
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